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. Abstract 

Oh' 

D I We analyze a class of mean-field (MF) lattice-fermion Hamiltoni- 

ans and construct the corresponding grand-canonical density operator 
for such system. New terms are introduced, which may be interpreted 
as local fugacities, molecular fields, etc. The presence of such terms 

Q ' is an unavoidable consequence of a consistent statistical description. 

Although in some cases (e.g. in the Hartree or the Hartree-Fock-type 

c/5 I approximation) the presented formalism is redundant, in general (e.g. 

for a renormalized t-J model) it leads to nontrivial modifications of 
the thermodynamic properties. 
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1 Introduction 

^ \ For most of realistic models of quantum many-particle systems an exact 

solution cannot be achieved. This is a typical situation condensed matter 

cn ■ or nuclear physics. As a consequence, a number of approximate methods 

have been developed. Among them, widely used are various types of the 
mean-field approaches. 

00 ! By mean-field (MF) description of the problem we understand making use 

^ \ of the Hamiltonian which depends on extra parameters, having the meaning 

of expectation values of well-defined operators. Values of those parameters 

^ . are not a 'priori known and are to be determined. The Hartree- Fock method 

is a good example of such approach. 

MF methods provide us with a valuable insight into the physical prop- 
erties of a system under the consideration. They are capable of describing 
many interesting phenomena, especially those involving phase transitions, 
such as that from normal metal to superconductor transition (e.g. in the 
form of the Bardeen-Cooper-Schrieffer theory), or those from a paramagnet 
to either antiferro- or ferro magnet. 
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In implementing such a procedure, it is the physics of the problem, not 
the mathematics, that tells us when the application of MF approach is jus- 
tifiable. Nonetheless, we may also be faced with the problem of considering 
the internal mathematical consistency of the method, as well as its underly- 
ing general principles. The aim of this paper is to clarify this issue to some 
extent, in a simple manner. 

The paper is organized as follows. Sec. [2]contains formulation of the prob- 
lem and invokes basic relevant notions from statistical mechanics. In Sec. [3], 
by means of the maximum entropy principle we construct the statistical- 
mechanical description of the system in the mean-field case. Explicitly, in 
I3.1l we derive the grand-canonical (GC) MF density operator, in 13.21 we com- 
ment on the MF thermodynamics and its consistence with the MF statistical 
mechanics, whereas in 13.31 we discuss the important case of the Hartree- 
Fock (HF) type MF Hamiltonians. Finally, in 13.41 the connection of present 
method with the variational principle based on the Bogoliubov inequality is 
elaborated. In Sec. H] two specific examples of mean-field Hamiltonians are 
provided, for which our approach leads to significant corrections. Section O 
contains a summary. 

2 Statement of the problem 

We assume that the Hilbert-Fock state space is of finite dimension D. The 
MF Hamiltonians are denoted as if, whereas those which are not of MF type, 
will be called exact and labelled with the subscript 'e'. 

As already mentioned in the Introduction, by mean-field approach we un- 
derstand the situation in which Hamiltonian depends on certain parameters, 
that are to be determined when solving the problem at hand. Such MF 
Hamiltonian is usually derived from the exact Hamiltonian Hf., by replacing 
some of the operators, {Ag} = Ai, . . . ,Am, appearing in He, by C-numbers 
{As} = Ai, A2, . . . , Am (= A)- Those numbers have a natural interpretation 
of expectation values of the corresponding operators {^s}. However, the 
mean-field Hamiltonian H may be not be related in an obvious manner to 
any exact Hamiltonian. 

We assume that the system is in thermal equilibrium with the parti- 
cle reservoir characterized by the temperature T = (3~^ and the chemical 
potential fi. The question then arises, how to determine the values of pa- 
rameters As appearing in H. If the Hamiltonian of the system were exact, 
i.e. not expectation-value dependent, then the treatment of such situation 
is straightforward, in the sense, that the grand canonical (GC) ensemble is 
well defined, and all information about the system is contained in the density 



operator [1, 2J defined as 

/3e = Z-iexp(-/3(^,-^iV)), (1) 

where Z^. = Tr[exp (— /3(ffe— /U^))] is the grand-canonical partition function. 
The expectation value of any observable A is then given by 

{A), = Tt[Ap,]. (2) 

Moreover, we assume, that in the mean-field case this basic definition still 
holds, i.e. 

{A) = Tr[ip(l)]. (3) 

Regardless of the form of MF density operator p, we should be able to obtain 
any expectation value through ([3]). In particular, the last definition allows 
us to determine the values of Ai,A2, ■ ■ . , Am appearing in p, by solving a 
system of M implicit equations for M variables (the so-called self-consistency 
equations), which take the form 

{As) = As = Tr[i,p(Ai, ^2, . . . , Am)], V,. (4) 

We may hope that this system should have at least one solution. 

One would expect, that MF grand-canonical density operator should be 
given by ([T]) with H^ replaced with H, so 

p = Z-^exp{-/3{H-^N)), Z = TT[exp{-l3{H-pN))]. (5) 

This is, however, not necessarily the case, as is discussed in detail below. 
In order to find the correct mean-field GC density operator for a given MF 
Hamiltonian H, we proceed analogously to the exact (non-MF) case and 
derive the form of such density operator from the basic principles of statistical 
mechanics. Thus, before concentrating in the remaining part of the paper 
on the mean-field case, we first would like to summarize briefiy the non-MF 
situation. Namely, for the exact Hamiltonian He, the GC density operator 
([T]) may be derived from the requirement that it maximizes the von Neumann 
entropy, S = — Trpelnpe, under the conditions that the expectation values 
of He and of the particle number operator A^, are fixed. This means, that we 
maximize the following functional (see [HE]) 

5e = Tr ( - Pe In Pe - PpeHe + (3ppeN - LOpe) ■ (6) 

The first two conditions are enforced with the help of the Lagrange multipliers 
(3 and p, whereas the parameter uj is introduced to ensure the normalization. 



Tr[/3e] = 1. Furthermore, p^ must be a constant of motion, so from the 
quantum Liouville equation [H [2] 

lh-p,{t) = [H,rPe{t)l (7) 

it follows, that [i5fe,Pe] = 0. If, additionally, [ffe,iV] = 0, which is usu- 
ally the case in non-relativistic quantum mechanics, we may simultaneously 
diagonalize each of the operators appearing in (El). 

Parenthetically, the maximization procedure is most conveniently carried 
out in the basis, in which pe is diagonal, i.e. pe = Ylii PjK)(*I- The resulting 
Pe is then of the form given by ([T]), [5i|. Density operators for the canonical 
and the micro-canonical ensembles may be obtained in a similar manner. 
The fact, that those density operators maximize the entropy functional ([6]) 
or analogous expressions, may be considered as their basic property [3]. 

3 Maximum entropy principle in the mean- 
field case 

3.1 Construction of the mean-field density operator 

We turn now to the construction of the proper form of MF density operator. 
We concentrate on the case, for which the entire system description is based 
solely on MF Hamiltonian. In the previous Section it was stated, that max- 
imum entropy condition is basic for canonical density operators. Moreover, 
the principle of maximum entropy is one of basic postulates of statistics, it 
ensures uniqueness and consistency of probability assignments, [3]. Conse- 
quently, it should be applied also to the particular case of mean-field descrip- 
tion. Now, the difference with the standard, i.e. non-MF situation above is 
that the Hamiltonian depends on M additional variables {^4^} = A, inter- 
preted as the expectation values of corresponding operators {^4^}. Those M 
expectation values, together with D diagonal matrix elements of p, are now 
the variational parameters, with respect to which the functional 

S = Tr( - plnp - m{A) + PppN - u{p - -^)), (8) 

is maximized. As earlier, —uj{p—l/D) terntJis added to ensure normalization 
of the density operator. Next, if those M + D parameters were treated as 



^Here taken in the form which shghtly differs from that of ^, but this technical 
alteration is inessential. 



independent, there is no certainty that {^4^} obtained in that manner coincide 
with those resulting form the self-consistency conditions (JH). In principle, one 
may use (jlj) to eliminate the Ag variables from H. Nonetheless, the resulting 
dependence of 5 on D diagonal elements of p is usually very complicated, 
and the MF density operator cannot be obtained in an explicit form. 

To resolve the situation, we use the Lagrange-multiplier method. Explic- 
itly, the M conditions (jl]) are enforced with the help of M Lagrange multi- 
pliers {As} = A. To account for that, the terms of the form As(Tr[pyls] — As) 
must be added to the functional ([H])- Under these circumstances, we have in 
total D+2M variables, {pi, . . . ,pd, Ai, . . . , Am, Ai, . . . , Xm}, treated as inde- 
pendent (apart from the normalization condition). In effect, the constrained 
entropy functional reads 



M 



Pxlnpx-p{pxH-fipxN-J2^sPxi^s-As))-uj{px--)^. (9) 

s=l 

The incorporation of the self-consistency constraints may be regarded as a 
redefinition of the MF Hamiltonian 

M 

H ^ Hx = H-Y,K{^s-As). (10) 

s=l 

The last step is the basic novel ingredient of our method and is discussed in 
detail in the remaining part of this paper. 

Few remarks are in place here. First, we are mainly interested in a sit- 
uation when each operator Ag is bilinear in creation and/or annihilation 
operatorqj, as then Hx is easily tractable. Second, some combination of the 
operators {Ag} appearing in (ITOl) may be proportional to the number opera- 
tor N. Nonetheless, the presence of the chemical potential p fixes the average 
number of particles N. 

Third, the operators H and iV in © have a different status than the op- 
erators {^s}. The expectation values of the former are assumed to be fixed 
a priori, enforced by the values of temperature and chemical potential. In 
contrast, the expectation values of the latter are obtained from a variational 
procedure, and are not a priori given. Consequently, the same distinction 
holds for parameters f3, p from one side and {A^} from the other, even though 
all of them play the role of Lagrange multipliers. Nonetheless, some of {As} 
may have a physical interpretation. Depending on the corresponding opera- 
tor As, they may be termed as local fugacities, molecular fields, etc. Note, 



^In principle, the presented formalism applies also to the general situation of the arbi- 
trary mean-field Hamiltonian. 



that for As = AJ, A^ is real, but for non-Hermitian Ag, it is a complex number 
in general. 

Apart from H and A^, there may be more observables with expectation 
values fixed a priori. In such a situation the terms {—^gAs), analogous to 
—fiN term, should be also added, to ensure that the quantity (Ag) is consis- 
tent with the prior information about the system. If we ignore some avail- 
able information about the physical system under consideration, we obtain 
a broader probability distribution. 

Fourth, the operators appearing in (ITOl) do not necessarily commute with 
each other, and thus cannot be, in general, simultaneously diagonalized. 
However, in analogy to the canonical case, we require that px thus obtained 
is stationary, i.e. dp\{t)/dt = 0, from which, using the quantum Liouville 
equation ([7]) for the MF case the condition [Hx,px] = follows. In particu- 
lar, in the case [A^, ^^a] 7^ (e.g. in the BCS theory of superconductivity), 
we can also redefine the Hamiltonian according to Hx — > Kx = Hx — pN, 
and then demand, that the corresponding grand Hamiltonian commutes with 
the density operator, [Kx,px] = 0. In result, diagonalization of Kx provides 
us again with the set of A and A- dependent eigenvectors {|«(v4, \))}fLi and 
corresponding eigenvalues of the density operator px- 

Fifth, the present formalism may be easily modified and applied to the 
case of canonical ensemble or discrete classical models (e.g. lattice gas, or 
spin systems). 

Finally, even if, strictly speaking, the T = situation cannot be inves- 
tigated directly within the framework of the present formalism, (must be 
analyzed separately), it is fully legitimate as well as convenient to analyze 
the situation in question as the T —>■ limit, e.g. for finite, but sufficiently 
low temperature. A more detailed analysis of this point will be presented 
elsewhere |S]. 

Let us go back to the problem of finding the maximum of the entropy sub- 
ject to the constraints. The functional ([9]), written in the common eigenbasis 
of Kx and px reads 

D M ^ 

Sx = ^l-qi\nqi-(3qi({H)ii-iJ,{N)ii-^Xs{{As)ii-As)) -uj{qi-—)\, 

(with {N)ii = {i\N\i) etc.) We assume, that the maximum of (TTTj) corre- 
sponds to stationary point. Then, the necessary conditions for an extremum 
of Sx with the normalization preserved, are 

dSx . A 






1-Z^'?.=0' (12) 



ds. 



oqj 



^^ s=l 



D 



(13) 



v.: l^-Z^E^lif-AO-O. (1*) 



and 



(9 An ^^" V9A 

1=1 



dS 



V^: TTT^ = /3X;^*((^-)--^"'))=0- (15) 



— * — * 

In ( TT^ and (IT^ , we ignored the possible explicit A- and A- dependence of the 
(normalized) eigenvectors \i). This is justified due to the Hellmann-Feynman 
theorem, |1]. Using (iT2l) and (iT3l) we obtain the explicit, basis-independent 
form of density operator p\, 

px = Z^'exp ( - f3{Hx - pN)), Zx = Tr[exp ( - (3{Hx - pN))]. (16) 



We see that (1161) differs from ([5]) by the presence of additional terms in the 
Hamiltonian, in accordance with (ITO|) . Additionally, Eqs. (IT^ and (IT5l) may 
be rewritten respectively as 

and 

Next, we define the generalized grand-canonical Landau functional 

^(1, A) = -p-' InZx = ~p-'Sx{pxiA, A), A, A), (19) 



with Zx given by ( lT6l) . That last relation may be easily verified by inserting 
Px given by (TT6!) into ([9]). Thus, the constrained minimum of J?^(A, A) corre- 
sponds to the constrained maximum of 5a(pa(^, A), A, A), here understood 
as a function of A and A only, once Eqs. ( IT2l) and ( IT3l) are solved. From now 
on we will use J^ instead of 5a. With the help of 1^, the Eqs. (^ and (HE]) 
are equivalent to 

Va^ = 0, Va^ = 0, (20) 



where Va = (g^; g^' ■ ■ ■ ' M~)' ^^^ analogously for Va- When deriving 
the last condition we made use of the identity (c.f. [2]) 



-e^W = / rf^e^^W ^dx) ]e^'-^^'-^^\ (21) 

ox Jq \ox / 

in combination with the linearity of the trace operation, as well as its invari- 
ance with respect to cyclic permutation of the operators. 

For each solution of (TT2 l) -( TT5l) or, equivalently, of ( l20l) . it should be verified 
separately if also the sufficient conditions for the existence of the constrained 
minimum of JF are fulfilled. In practice, this will be most conveniently carried 
out through numerical analysis. 

We may expect that one of the solutions of (^01) . labeled (Aq, Aq), yields 
such desired global minimum of JF subject to constraints, and this defines 
the equilibrium state. If more than one solution with the minimal value of 
JF exist, we select one of them, breaking the remaining symmetry. 

However, JF is well defined also for the values of A which differ from 
Aq, and may possess physical interpretation provided that the consistency 
conditions ( ITSl) are met. Namely, J^{A, \{A)) may be interpreted as a Landau 
functional corresponding to the given MF model. The explicit analytical form 
of jF(yl, A(A)) or its series expansion in the powers of order parameter(s) 
may be difficult to obtain, and then again we have to resort to the numerical 
analysis. Nonetheless, the concepts of Landau theory of phase transitions 
may be utilized. E.g., the quantity A/""^ exp ( — /5jF(y4, A(y4))), where M is 
normalization constant, may be interpreted as the probability distribution 
of A. Also, plots of T{A, \{A)) may be used for qualitative analysis of the 
character of the emerging phase transition. 

3.2 Mean-field equilibrium thermodynamics 

The equilibrium description is justified, provided that at the point A = 
Aq the function A/'~^exp( — pj^{A, \{A))) possesses a sharp peak, and 
consequently, the fluctuations of order parameter can be neglected. This 
condition is expected to hold for a wide class of MF Hamiltonians deflned 
on lattice, in the A ^ oo limit, where A is the number of lattice sites. 
Obviously, for such models, this condition is equivalent to the thermodynamic 
limit V — ^ cxo, N ^ oo with finite N/V = n. Then we can define the 
thermodynamic grand potential 

fi(T,V,/i)=^(r,V,/i;lo(r,y,/x),Ao(T,V,/i)). (22) 

With this definition, the MF thermodynamics can be constructed, but we 
still have to demonstrate its consistency with MF statistical mechanical de- 
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scription based on flTUl) . cf. [6j. Such consistency is certainly present for an 
exact description, making use of non-mean-field Hamiltonians [H El El E] ; in 
that case we have for example 

^"-feJv.- P=-(i^)r.- ""--{I^ry- ^''^ 

One may ask if in the MF case the r.h.s. of (|23l) correspond to their statistical- 
mechanical definitions, as the equilibrium MF Hamiltonian Hx obviously 
depends on temperature and chemical potential through its dependence on 
the equilibrium values of mean fields and Lagrange multipliers^, (y4o,Ao). 
It can be shown, that fl^ indeed holds, due to the underlying variational 
principle. Explicitly, 



idf)v,, ^ (3%{(3-HnZ,)=(3Z^'dpZ,-\nZ, 



A 



= -pZ~'Ti[{pdpKx + kx)e-^^']-\nZx 

= -(3^{d^Kx)x-(3{Kx)x-\nZx 

= -(3^d^Kx)x + {\npx)x = ~/3'{df,Kx)x~S, (24) 

where Kx = Hx — iiN and Zx are evaluated for A = Aq, A = Aq. Also, 
hiZx{AQ, Aq) = — /3f^, and {Kx)x = U — fiN. In effect, we obtain the desired 
result fl2H|) . provided that the term {dj3Kx)x vanishes. Indeed, in equilibrium 
we have 

/a T>\ /^^^ fdkxdAso , dkxdXso 



s=l 

M 



s=l 



dA, dp dXs dp 

dAso / dHx \ 
X dp \ dXs Ix dp 



^ V\ dA, /x dB \ dX, /x dl 



(25) 



due to relations (1171) and (1181) . The same reasoning applies to the derivatives 
of Vt with respect to p and V . In the latter case, Hamiltonian depends 
explicitly on V (e.g. through the volume dependence of its parameters). 
In effect, we recover proper expressions for average number of particles N 
and pressure p. Other thermodynamic potentials may be constructed in the 
standard manner, e.g. F = Vt + fiN. 

Note that as the P- and /i- dependences of equilibrium values of mean- 
fields Aq obtained through self-consistent Eqs. (jl]) are nontrivial in general. 



We exclude the possibility of explicit dependence of H\ on f3 = \/T or /i. 



and also because fITTj) with A = usually does not hold in such case, there is 
no reason for {dpK) (as well as for analogous expressions for the derivatives 
with respect to yU, V etc.) to vanish. This feature results in the lack of 
consistency between MF thermodynamics and MF statistical mechanics even 
on the level of the first-order derivatives of the corresponding thermodynamic 
potential and constitutes an obvious and serious drawback of non- variational 
self-consistent approaclu. The problem of such consistency is also analyzed 
in Ref. [6], however, from a slightly different perspective. The analysis of 
the second-order derivatives of Vt is slightly more complicated and is carried 
out in [5]. 

3.3 Redundancy conditions 

The natural question is, when the presence of Lagrange multipliers in the 
Hamiltonian ( ITO|) is unnecessary. We will not provide the general answer, but 
instead present an important example. Suppose we have the exact Hamilto- 
nian of the following form 

H, = Hoe+Y, y-iO-O'y^ (26) 

{k,7} 

where each O^ is assumed to be bilinear in the creation and/or annihilation 
operators, and the summation is taken over all pairs {/t, 7} of multi indices. If 
the interaction term in (l26l) is decoupled (i.e. replaced by its MF counterpart) 
according tqfl 

6^6^ -^ A,At + A,At - A,At, (27) 

with As, At given by (jl|) for the appropriate density matrix, then the standard 
Hartree or the Hartree-Fock (HF) type mean-field Hamiltonian is obtained. 
With the constraint terms it reads 

^A = HQe + ^Vst{AsAt + AsAt- AsAt) -^\s{As- As) 

{s,t} s 

= i^-5^A,(4-^). (28) 



"^A notable exception, the case of pure Hartree or Hartree-Fock type of MF Hamiltoni- 
ans, is discussed in the next paragraph. 

^Usually, the O^Oj term may be decoupled in more then one way, but this does not 
change the arguments presented. 
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In the above, i^oe does not depend on A, and Kjt are the combinations of 
matrix elements Vk^. Then, one can find that 

^ s 

which together with flT7|) imphes that A^ = 0, for any t. Also, for Hamiltonian 
of such form, the solution of M self-consistent equations (jl]) coincide with 
those of M equations Va^ = of unconstrained variational approach. A 
related discussion of this problem may be found in [7J. 

Summarizing, for any MF Hamiltonian of the Hartree-Fock type, applica- 
tion of our method is not necessary. However, there are important examples 
of the MF Hamiltonians not being of HF form, in Section H] we present two 
such examples. 

3.4 Relation of the present method to the variational 
principle of Bogoliubov and Feynman 

In the previous Sections we considered the situation, when the description 
of the system has been based solely on the MF Hamiltonian H. In other 
words, we made no use of any exact Hamiltonian He present in a problem at 
hand, from which H may be derived and for which it may play the role of a 
simplified counterpart. 

However, frequently the MF density operators are used as a trial varia- 
tional states in order to determine the bounds on the grand potential Qg or 
free energy Fg of the system described by He- Such bounds are provided by 
the well-known inequality due to Bogoliubov and Feynman [8], [9] 

fie < {Ke - A^o)o + fi = Tr[po(A:e - ATo)] + fi. (30) 

In the above, Hq is a trial Hamiltonian, Cl = —P~^ lnTr[exp (^—P{Hq—^qN)^], 
and Qe = -/3"MnTr[exp ( - /3{He - fieN))]. Also, K^ = H^ - fieN and 
analogously for Kq, with /ig (/xo) being respective chemical potentials. There 
are several ways to prove (l3Up . probably the shortest of them is that of making 
use of the Klein inequality, [TU] 

Tr[plnp] > Tr[plna], (31) 

which holds for any normalized density operators p and a. Inserting 

p = Po = Zq^ exp{-(3{Ho - fioN)), a = pe = Z^^ exp(-/3(^e - P'eN)), 

(32) 
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(with Zo = Tr[exp(-/3(ifo - f^oN))] and Z^ = Tr[exp(-/3(ife - fJ^N))]) into 
flHTj) . we obtain flHUj) . If one wants to work with the chemical potential as an 
independent variable, then /Xg = /io = /^ and ( !30|) reduces to 



fie < (i^e - Ho)o + n = TT[po{He ' ^o)] + ^- (33) 

However, usually it is more convenient to have particle number A^ as an 
independent variable. In such a case /io and (unknown) /ig are not equal in 
general, as they are determined from different conditions, i.e. (iV)o = N and 
(iV)e = A^, respectively. However, flHUj) may be then given the form 

Fe = fie + fieN < {H, - Ho)o + fi + fioN. (34) 

Bogoliubov inequality is the basis of the variational principle, the best trial 
state po is the one minimizing r.h.s. of (IH^ . or (IMIl . depending on the choice 
of independent variables, /i or A^, respectively. 

If one wants to use the Bogoliubov inequality in order to determine upper 
bounds for fig or F^ of some Hg, but with the self-consistency of the approach 
being preserved, then the density operators ([5]) are not good candidates for 
the trial states. Clearly, they have no free variational parameterq^ if only 
the values of mean- fields, A = Aq, are obtained in a self-consistent fashion 
through (jl]). On the other hand, if Aq were obtained from unwary minimiza- 
tion of the r.h.s. of ( l33l) . ( l34l) with p given by ([5]), then the self-consistency of 
the approach would be violated in general (again with the notable exception 
of the Hartree-Fock MF Hamiltonians) . 

In contrast, consider the trial state chosen in the form 0161) . with the trial 
Hamiltonian Hq = H\, yet unrelated to the He in question. Then the r.h.s. 
of (l33l) reads 

Ti[pxiHe - Hx)] + ^(i. A) = BiA, A), (35) 

whereas that of flH^ . respectively (with po = p) 

Tr[pA(^e - Hy)] + ^(1, \)+pN = B{A, X) + pN = i3^(l. A). (36) 

Depending on the choice of independent variable, the necessarju 2M condi- 
tions for r.h.s of the Bogoliubov inequality, i.e B or Bn, to have a minimum 
subject to constrains of self-consistency acquire slightly different form. For 
(135|) we have 

VaS(A, A) = 0, Va^(A A) = 0. (37) 



^We neglect the situation, in which there appear some extra variational parameters, 
not being expectation values of operators present in a problem, and thus not requiring the 
corresponding Lagrange multipliers. 

^Similarly as previously, here we assume, that the desired minimum corresponds to a 
stationary point. 
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On the other hand, for flHbl) we have 

VUBn{A, A) = 0, Va(-F(A A) + fiN) = 0, ^iHA A) + fiN) = fi. 

(38) 
In the above, V'^ labels the gradient taken with respect to all mean-fields 
except Ai = N. Also 

V^Bn{A, A) = V^B{A, A), VxiHA A) + fiN) = VxHA A). (39) 

If B{A, A) 7^ J^{A, A), the solution of ( 1371) or ( l38l) leads to different values 
of Aq and Ao, than those obtained from (1201) . Then, p^ is no longer a true 
grand-canonical MF density operator for Hx- The above remarks hold for 
other inequalities (i.e. various generalizations of (1301) ). and the corresponding 
variational principles built on them, [9]. One may use the px (1161) as useful 
ansatze in such cases as well. Nonetheless, the entropy functional ([9]) with the 
corresponding variational principle plays a special role, and it determines the 
form (ITUj) of GC MF density operator. Also, the self-consistency conditions 
may be conveniently formulated with the help of J^{A, A), as in (1H7|) and (IH^ . 
However, clearly, for any pair of He, Hx, for which the condition 

{H,-Hx)x = (40) 

holds, so consequently also Q{A,X) = J^{A,\) = B{A,\), the Bogoliubov - 
Feynman and Maximum Entropy variational principles coincide. In such case 
the optimal, from the point of view of maximum entropy inference, i.e. the 
most uninformative MF density operator px provides us also with the upper 
bound for exact grand potential Qe or free energy F^ of Hg. For (1371) this 
is obvious. For (|38l) . this may be shown using dJ-'/dfi = —N and replacing 
B ^ J-", then we obtain the identity 



Consequently, dJ-'/dN = for the solutions of Eqs. ( l38l) . From Eqs. ( l39l) 
and (HTl) . it follows that any solution of (l38l) . i.e. {Aq, Aq, p) is also a solution 
of (IHTj) . if A^ and p related by the condition (iV) = A^ = —dT /dp. Also, 
conversly, Eqs. (IHUj) . (HTl) imply that the solutions of (IH71) are also solutions 
of (!38l) . thus both schemes are equivalent. 

4 Two nontrivial examples 

We now discuss briefly two MF models, for which the proposed method leads 
to nontrivial corrections. To show this, we make use of Eq. ( IT71) . and argue, 
that the condition A = cannot be, in general, satisfled. 
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4.1 Renormalized t-J model 

The first example we analyze is the MF Hamiltonian of the so-called renor- 
malized t-J model ^ilj-^iSj. It originates from the standard t-J model [111, ^^~ 
lieved to describe correctly the essential physics of the cuprate high-temperature 
superconductors, that is expressed by the Hamiltonian 

i,j(i),a (ij) i,cr 

Here Cjo- = (1 — ni_^)ci„ and z/j = ^^ = (y^^Cifj are the operators with double 
occupancies on site i projected out. The corresponding renormalized MF 
Hamiltonian may be taken in the form |13i il5j 



^ = XI [{tijSlflaCja + H.C.) - -Jgij{Xjic\aCj^ + H.C. - \x^j\^) 

+ -^ j4(A,,4c1_, + H.C. - I A,,f )) - /i J] cLq., (43) 

with rii = J2a{4aCia), Xij = {4a(^ja) , aud Aij = {Ci^aCja) = {Cj-aCia) ■ In 

this form the projections have been abolished at the price of introducing 
the expectation- value dependent renormalization factors g^j and g:/j resulting 
from the Gutzwiller approximation [I8]|. At first, we consider the renormal- 
ization factors depending solely on Xi = 1 — rii, i.e. 



The task is to determine the values of the mean-field parameters Xijy ^ij and 
Xi. This could be achieved with the help of self-consistent Eqs. (jlj) for H 
given by (|43]) (cf. [Il]- [IS]). 

On the other hand, we may proceed in the framework of our method, and 
redefine the Hamiltonian (I43p according to 

Hx = H - ^Ar(5^cLQ.-n.)-5^(Af,.(cU.-x^,) + H.c.) 

i a {ij)a 

Assuming no additional symmetries, using d/dui = —d/dxi, converting the 
sums over bonds to that over sites in (HHj) according to 2 J^u ,) = J2i jU) ^^^ 
finally taking A^, = rii in Eqs. ( TTTll we obtain 
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2^" = (E^^.%(4^-- + H-c.) 



j{i),'^ 






Y, (««^Rex„ - l^^d*/ + |A./)). (46) 



Below it will be shown that in this particular case, A^- = A^ = 0, for each 
bond (ij). Consequently, if for all i, we put A" = 0, then the density opera- 
tor, and all averages coincide with those of the standard self-consistent MF 
treatment. However, in such situation equations (HUj) and the self-consistent 
equations (jll) for H given by ( H3l) cannot in general be simultaneous satisfied. 
For example, if one sets J = and Xi G (0, 1), then from (H6l) it follows that 



-.t 



E(^<=.U. + H.c.))=0. (47) 

i(«)o" 

This is a senseless result, as in such case our Hamiltonian reduces essentially 
to that of free fermions, at least for the case of homogeneous solution, Xi = x. 
Similarly, for tij = 0, one obtains |xjjP = 0, |Ajjp = for each bond ij, 
provided that each tij is of the same sign, as dgjj/dxi > 0, for any Xi and 
(ij). This is in contradiction with the numerical results of Raczkowski et ai, 
|14j . In other words, in this particular case, we can infer about the predictions 
of our method using the previous, non-variational results. Consequently, it 
follows that A" 7^ 0, for every i. 

From the point of view of our findings, a related situation can be found 
in Ref . [16] . Namely, the inhomogeneous solutions of the model (1431) (with 
slightly different definitions of Xij ^ind Ajj), and with the renormalization 
factors given by (jH]) are analyzed. However, the additional terms of the 
form X]i^«(X]cr'^Lcjo- + Xj — 1), are included, where e^ is the local fugacity. 
The quantity e^ plays a role analogous to A^ in our method. Inclusion of that 
additional constraint allows to treat each Xi as a variational parameter, and 
leads to the equations analogous to (HUjl . This should be clear in the light of 
the discussion in the preceding Sections. 



15 



Next, we pose the question of redundancy of the Lagrange multiphers Xfj 
and A^. We may expect, according to discussion carried out in Section [X^ 
that these quantities vanish, as the Hamiltonian ( H3|) with g^j and g^j given 
by (jSl) is of the Hartree-Fock type with respect to both Xij and Ajj. This is 
indeed the case, we may use once more Eq. ( TT71) with A^j = Xij or A^ = Ajj, 
respectively, to obtain 

^g = (i^E4(44-.-^^)),-0- (49) 

Obviously, those self-consistency equations must be fulfilled both in the pre- 
vious treatments (A = 0) as well as in the present method. 
Next, we consider the renormalization factors of the form 



t _ I 4:XiXj[l Xi)[l Xj) , . 

,J ^ A{1 - Xi){l - Xj) 

^'^ (l-xf)(l-x2) + 8x,x,(|A,,f-|x.,f) + 16(|A,,|4 + |x.,f)- ^ ' 



In such more general situation, the equation fHbjl is still valid, and Eqs. 
and f H9l) acquire the form, respectively 



^^ = -E(^^.#(4^- + H-cO-I^M(x.4.c,. + Hx.-|^ 



r 

in/ 

, ...^,j ' 2" dxi 



+ -7^4(4^- - X^) - \j^{/^.,c]A-. + H.C. - |A, '^^ 

-4t,,|^Rex., + ^j|^(|A.,f + |x.,f), (52) 



-^S = E(l^M(x^-^4^^'^ + H-c.-|xd') + !^4(4-.4-^*.) 



^/'A..^t J I U „ . IA..|2^ 



£j_^fA-c^ c^ +Hc-IA- 
4 9A- -^ ■^ 

^M(|A,f + |x,f)- (53) 



2 9A,j 
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Again, if we put A = in order to reduce our approach to the standard 
treatment, that Eqs. fl3^ and fl3^ cannot be satisfied. For Eq. (1321) this can 
be easily shown for the case tij = ~t by adding it to its complex conjugate 
part. Then, at least for the class of solutions with for Rexij > 0, one can 
convince oneself by explicitly computing respective derivatives of g\j and gf^, 
that the resulting expression should have nonzero value, for a wide range 
of values Xi. For Eq. (13H11 or (1321) in the case of tij = 0, no additional 
explanation is required. 

Similar arguments apply if, instead of A = 0, a weaker conditions, A^ = 
or A^ = 0, are imposed. In such situation we can no longer make use of 
the numerical or analytical results of the purely self-consistent approach. 
Nonetheless, it is unlikely that such conditions do not lead to contradictions. 
E.g. from Eqs. (I32|) with t^j = or from (13^ in a general case, we infer, 

that Xij = ^ij = 0. 

In that slightly hand-waving manner we tried to argue, that for the Hamil- 
tonian (H3l) our variational approach does not reduce to that based solely on 
the self-consistency equations. 

This model has been recently analyzed and solved numerically in a sepa- 
rate publication [19]. 



4.2 Electrons with spin-dependent effective masses 



In Ref. [20], an approximate one-particle description of a strongly-correlated 
electron system was introduced, based on the observation, that the effective 
masses of quasiparticles may be spin-dependent. The MF Hamiltonian of 
this model depends on Ai = N and A2 = M, where 

^ = ^1 = E E ^i^ck., M = A2 = Y,Y1 ^^Ic^ck.. (54) 

k cr k IT 

With additional constraint terms, the Hamiltonian reads 

ffSDM^ (5^5^($.(iV,M)ek.-a/i,-/xK.)-Ai(iV-iV)-A2(M-M), 

k (T 

(55) 
where {—"^^crhanka) is the Zeeman term, and ha is the reduced applied 
magnetic field. The present case is, in fact, quite similar to the previous one, 
the band-narrowing factors ^criN,M) are also derived using the Gutzwiller 
approach. Then analogously, by making use of Eq. (fT7|) . either with A^ = A^ 
or A^ = M, we obtain respectively 

/ g$,(iV,M) . \ 
Ai = -(2.2. g^ ek.nk.)^, (56) 

k cr 

17 



EE 



0^ ek^^ka).. (57) 



Regardless of the detailed analytical form of ^a{.N, M), we multiply (IFr|) by 

and add it to ( l56l) . The resulting r.h.s is then proportional to (i?fcm,t);^ = 
( Sk^kt^kT);^) the average band energy of particles with 2;-spin component 
a =t. If we put Ai = A2 = 0, all the averages become identical to those of 
the standard treatment. Consequently, for (13U1) and (IFH) to hold in such a 
case we have to require that (-Efcin.t) = O5 fo^ ^iny admissible value of N, M, 
and functional dependence of $o-(A^, M) on its variables, which is a trivial 
and inconsistent result. Obviously, the same holds for o" =],. 

The numerical studies of this model and its extensions are carried out 
recently. 



5 Summary and conclusions 

In this paper we have presented the method, based on the maximum en- 
tropy principle, of constructing the grand-canonical density operator for a 
wide class of mean-field Hamiltonians. We have shown, that such a density 
operator is not, in general, obtained by replacing the exact Hamiltonian by 
its mean-field counterpart in the grand canonical density operator. Usually, 
some extra terms, which may be interpreted as a kind of source terms or 
molecular fields, must be added to the mean-field Hamiltonian. This mod- 
ification also ensures consistency of basic thermodynamic relations, which 
is not guaranteed for self-consistent methods which are not based on the 
variational principle. 

It is also shown, that although for the Hartree-Fock type of mean-field 
Hamiltonians application of our method is not necessary, there are also some 
important examples, for which it does not reduce to the standard approach. 
The renormalized mean-field t-J Hamiltonian is such case. 

If there is M expectation values appearing in a mean-field Hamiltonian 
that are to be determined, the method presented here results in a system of 
2M equations. This is in contrast to the standard treatment, where only M 
such equations appear. 

The problem of the application of maximum entropy inference in the 
mean-field case was studied also in a series of the papers [22]. However, the 
presented point of view and the results obtained differ essentially from ours. 
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Also, the problem partly related to that discussed in the present paper, 
namely the construction of variational principles suited for optimization of 
physical quantities of interest, is examined in an interesting article [23]. This 
paper contains a detailed analysis of the problem, as well as covers a large area 
of the subject (for example the authors discuss time- dependent formalism). 
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